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LONG-TERM  GOALS 


The  goal  is  to  bring  to  fruition  R&D  that  has  been  performed  for  ONR  since  FY05.  The  work  has 
involved  the  development  of  a  state-of-the-art,  high-fidelity,  high-speed,  3-D  finite-element  (FE) 
broadband  computer  simulation  capability  for  modeling  the  scattering  of  sonar  signals  by  undersea 
objects  located  in  or  near  the  seabed  in  littoral  environments  with  smooth  or  rippled  water/sediment 
interfaces  (Fig.  1).  Many  of  the  computational  techniques  developed  in  this  effort  will  be  relevant  to 
future  ONR  programs  involving  different  applications. 
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Figure  1.  A  representative  shallow-water  acoustic  scattering  scenario. 


OBJECTIVES 

The  primary  objective  is  to  provide  the  U.S.  Navy  with  a  high-fidelity,  high-speed  computer 
simulation  tool  for  computing  and  displaying  the  frequency  and  aspect  angle  response  (a.k.a.  “acoustic 
color”)  of  complicated,  realistic  objects  in  littoral  environments.  A  secondary  objective  is  to  maintain  a 
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continuing  R&D  effort  in  order  to  significantly  increase  the  computational  efficiency  of  the  software, 
thereby  expanding  the  ability  to  model  more  complicated  objects  at  higher  frequencies. 


APPROACH 


>  High  fidelity 


Physics:  Target  scattering  models  employ  fully  3-D  physics  throughout  the  computational 
domain;  i.e.,  the  fully  3-D  wave  equations  are  used  throughout  both  solids  and  fluids; 
Fluids:  linear  acoustics  (Helmholtz  eq.):  inviscid,  volume  dissipation,  inhomogeneous 
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Solids:  linear  elasticity  (elastodynamics  eq.):  viscoelastic,  anisotropic,  inhomogeneous 

-V.(cV«)-  (o2psu  =  fs 


(2) 


Equation  (2)  is  used  throughout  all  structural  parts  of  the  target,  including  even  thin  plates  and 
shells,  small  fillets  and  chamfers,  etc.  This  is  essential  for  high  fidelity  since  the  physics  is 
strongly  3-D  near  discontinuities  such  as  shell  intersections.  Local  evanescent  elastic  waves  must 
be  modeled  accurately,  since  their  phases  and  amplitudes  detennine  the  propagating  elastic 
waves  in  the  outer  casings,  which  in  turn  determine  the  scattered  acoustic  field  launched  from  the 
casing. 

•  Mathematics:  Since  the  infinite  family  of  polynomial  basis  functions  within  each  finite  element 
spans  a  complete  space,  then  the  FE  solution  within  each  element  can  be  arbitrarily  close  to  the 
exact  solution  of  the  physical  model. 


In  summary,  all  physical  phenomena  describable  by  eqs.  (1)  and  (2)  are  included  in  the  FE 
solutions. 


>  High  speed 

Acoustic  color  is  target  strength  (TS)  as  a  function  of  both  frequency  and  aspect  angle.  To  generate  an 
acoustic  color  plot  using  computer  simulation  requires  extraordinarily  high  computational  efficiency, 
since  it  involves  a  double  sweep;  (1)  over  a  broad  frequency  band,  typically  several  octaves,  which 
usually  entails  sampling  at  several  hundred  frequencies,  and  (2),  for  each  frequency,  a  sweep  over 
many  aspect  angles,  again  typically  several  hundred.  The  double  sweep  therefore  requires  O(105)  FE 
models.  In  addition,  by  decomposing  the  incident  field  (Fig.  3  under  Work  Completed)  into  different 
types  of  symmetry  components,  each  model  can  often  be  decomposed  into  2,  4  or  8  smaller  (hence 
much  faster)  sub-models.  The  net  result,  therefore,  can  sometimes  involve  O(106)  FE  analyses  to 
produce  a  single  acoustic  color  plot.  For  practicality,  the  goal  of  this  work  is  to  compute  an  acoustic 
color  plot  in  about  one  day,  which  means  Q(106)  FE  models  per  day.  This  goal  is  close  to  being 
accomplished,  by  a  combination  of  software  and  hardware  development,  as  described  below  under 
Work  Completed. 
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>  High  reliability 

The  control  of  errors,  via  verification  and  validation  (“V&V”),  plays  a  key  role  in  this  project;  it  is 
essential  to  the  goal  of  providing  reliable  simulations  of  the  real  world.  Verification  (Fig.  2)  refers  to 
the  control  of  errors  in  the  finite-element  solution  vis-a-vis  the  exact  solution  of  the  idealized 
mathematical  model,  i.e.,  computational/mathematical  errors.  Validation  refers  to  the  control  of  errors 
in  the  exact  solution  of  the  idealized  mathematical  model  vis-a-vis  the  real  world,  i.e.,  errors  in  the 
physical  assumptions  and  data  in  the  idealized  model.  V&V  is  an  ongoing,  never-ending  exercise  as 
the  simulation  capabilities  mature. 
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Figure  2.  Pictorial  explanation  of  verification  and  validation. 


References  [1]  and  [2]  illustrate  several  verifications.  Figure  5  under  Results  below  shows  a  significant 
validation  milestone:  experimental  and  simulated  acoustic  color  plots  exhibiting  a  striking  similarity. 
V&V  is  an  ongoing,  never-ending  exercise  as  the  simulation  capabilities  mature. 

>  Key  individuals  participating  in  this  work 


•  Leslie  H.  Wigdor:  Syslink  Consulting,  LLC.,  21  Rock  Street,  Cold  Spring,  NY  10516-2911. 
Systems  programming;  integration  of  FE  and  non-FE  software;  user-interfaces 

•  Dr.  Gary  S.  Sammelmann:  Code  HS1 1,  Naval  Surface  Warfare  Center,  1 10  Vernon  Ave., 
Panama  City,  FL  32407-7001.  Development  of  non-FE  theory  and  algorithms;  analytical 
solutions  to  canonical  problems  for  verification. 

WORK  COMPLETED 


>  Finite-element  R&D 
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As  reported  in  [1]  and  [2],  several  modeling  techniques  were  developed  during  FY06  and  FY07: 

•  New  tri-axial  ellipsoidal  radiation  boundary  conditions, 

•  Group-theoretic  decomposition  of  the  incident  field, 

•  Parametric  sweeping  over  both  frequency  and  aspect  angle, 

•  Evaluation  of  scattered  field  exterior  to  the  FE  mesh  using  Kirchoff-Helmholtz  integral. 

The  motivation  for  developing  all  these  techniques  has  been  to  increase  computational  efficiency.  Each 
technique  has  increased  computational  efficiency  by  one  to  two  orders  of  magnitude.  The  resulting 
half-a-dozen  orders-of-magnitude  increase  in  computational  efficiency  has  resulted  in  meeting  the 
high-speed  requirement  of  O(106)  FE  models  per  day. 

During  FY08  these  techniques  were  further  refined.  However,  the  major  effort  was  integrating  them 
into  an  acoustic  color  system,  which  is  described  next. 

>  Software:  automated  acoustic  color  system 

Most  of  the  work  in  FY08  has  focused  on  bringing  together  all  the  above-described  technologies  into  a 
user-friendly  GUI-driven  software  system  capable  of  automatically  running  O(106)  FE  models, 
managing  computer  resources  and  very  large  amounts  of  different  types  of  data,  and  combining  the 
output  from  all  the  models  into  an  acoustic  color  plot,  as  well  as  a  variety  of  other  plots,  e.g.,  slices  of 
the  acoustic  color  representing  frequency  response  or  angular  response.  The  principal  functions  are 
listed  here,  followed  by  an  illustrative  plot  (Fig.  3)  of  a  scattering  problem  that  is  modeled  by  only  one 
quadrant,  due  to  symmetries,  and  a  sketch  of  the  surface  on  which  Helmholtz  integration  is  performed. 

•  Sweeping  over  frequency, 

•  Sweeping  over  aspect, 

•  Management  of  sub-models,  due  to  symmetries  in  target  and  environment; 
the  incident  field  can  be  nonsymmetric, 

•  Dynamic  frequency-dependent  geometry  and  meshing  to  maintain  uniform 
discretization  error  with  frequency  -  an  important  feature  for  signal  processing  the 
acoustic  color  data.  (See  illutration  and  explanation  under  Results  below.) 

•  Computation  of  Helmholtz  integral  for  TS: 

-  evaluating  near-field  p  and  dp/dn  on  Helmholtz  surface,  which  includes  “virtual” 
surfaces  that  don’t  exist, 

-  evaluating  G  and  dG/dn  for  littoral  environment, 

.  GUI: 

-  Input:  parametric  ranges,  increments,  etc., 

-  Output:  Acoustic  color  data,  in  various  plots. 
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Figure  3.  A  typical  model:  one  quadrant  of  a  full  model  of 
scattering  from  a  buried  spherical  shell.  The  dashed  contour 
indicates  the  edges  of  a  closed  surface  for  Helmholtz  integration. 

Figure  4  illustrates  the  user-friendly  GUI  for  controlling  a  complete  acoustic  color  simulation. 


Figure  4.  GUI  that  controls  production  of  acoustic  color  plots  and  other  output 
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>  Hardware:  building  a  distributed-processing  rack-computer  system 

As  demonstrated  under  Results  below,  an  acoustic  color  simulation  can  typically  take  a  few  processor- 
days  for  “practical  accuracy”  (~  1  dB  discretization  errors)  or  as  many  as  several  tens  of  processor- 
days  for  “high  accuracy”  (~0.1  dB  errors).  Distributed  processing  is  therefore  necessary  to  achieve  the 
high-speed  requirement  of  approximately  one  day  to  compute  an  acoustic  color  data  set.  An  important 
task  during  FY08  has  been  the  design,  purchase  and  assembly  of  a  scalable-architecture,  multi-blade, 
parallel-processing  rack  computer.  This  is  almost  finished,  including  the  installation  and  successful 
testing  of  the  acoustic  color  system  software.  When  completely  finished  (in  1QFY09),  the  system  will 
have  20  blades,  each  with  two  dual-core  processors  (~  3x  speed-up  vs  one  core).  Ten  blades  will  have 
32  GB  RAM  (for  high  frequencies)  and  ten  will  have  16  GB  RAM  (for  medium  and  low  frequencies). 
This  combination  of  blades,  RAM  and  multi-core  processors  will  reduce  several  tens  of  processor  days 
to  about  one  day. 

In  parallel  with  the  hardware  build  has  been  the  task  of  modifying  the  acoustic  color  system  to  manage 
computer  resources  in  a  distributed-processing  environment,  with  concomitant  modifications  to  the 
GUI. 

RESULTS 

An  important  milestone  was  achieved  with  the  following  experimental  validation.  The  acoustic  color 
system  was  used  to  model  the  free-field  insonification  of  a  5-to-l  aspect  ratio  aluminum  cylinder 
(length  7.5”;  diameter  1.5”)  from  50  to  120  kHz  ( ka  from  4.2  to  9.7;  kL  from  42  to  97).  The  results  are 
compared  here  (Fig.  5)  with  the  data  from  a  tank  experiment  performed  at  Wash.  State  Univ.  [3],  This 
was  a  “practical  accuracy”  simulation,  using  a  relatively  coarse  sampling  in  frequency  and  aspect 
angle,  which  required  “only”  102,648  FE  analyses.  A  “high  accuracy”  simulation,  using  finer 
sampling,  would  require  O(106)  FE  analyses  and  take  twenty  to  thirty  processor  days,  i.e.,  about  one 
day  on  the  above-described  distributed-processing  rack  computer. 
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•  A f=  1  kHz  ;  A<?>=  1° 

•  Discretization  error  ~  1  dB 

•  Average  mesh  size  =  83,000  dot 
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Figure  5.  FE  model  vis-a-vis  experimental  data 


This  validation  problem  illustrates  the  use  of  the  automated  dynamic  change  to  both  the  outer  fluid 
boundary  radius  and  the  mesh  density  throughout  the  model.  As  frequency  increases:  (1)  the  outer 
fluid  boundary  automatically  moves  in  to  maintain  a  constant-wavelength  separation  from  the  target 
and  (2)  the  mesh  is  automatically  refined  to  maintain  a  constant  number  of  elements  per  wavelength. 
Figure  6  shows  the  FE  models  at  three  of  the  141  frequencies  run.  The  meshes  are  on  one  of  the  faces 
of  the  3-D  model  in  the  upper  right  of  the  figure,  which  is  one  octant  of  the  full  physical  model. 
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Physical  model:  FE  model: 

Symmetries  =>  only  1  octant  is 
oo  necessary,  for  all  aspect  angles: 


Figure  6.  Automated  change  of  geometry  and  mesh 
with  frequency,  showing  3  of  the  141  different  FE  models 


IMPACT/APPLICATIONS 

This  high-speed,  high-fidelity  3-D  target  scattering  simulation  capability  will  be  an  important  asset  for 
SWAMSI  (Shallow- Water  Autonomous  Mine  Sensing  Initiative),  BMC  (Broadband  Mine 
Classification),  LFBB  (Low  Frequency  BroadBand)  and  the  SERDP  UXO  (Strategic  Environmental 
Research  and  Development  Program,  UneXploded  Ordnance)  Thrust  Area,  since  it  is  now  possible  to 
generate  acoustic  color  data  sets  for  many  objects  and  environmental  scenarios.  The  scientific  impact 
will  be  all  the  new  modeling  techniques  developed,  which  should  be  useful  in  other  DoD  programs  as 
well  as  in  the  general  structural  acoustics  community. 

In  FY04  and  FY05  Dr.  Burnett  researched  the  commercial  FE  software  and  decided  to  use  the  code 
Comsol  Multiphysics  (then  called  Femlab)  for  modeling  acoustic  scattering  from  objects  in  shallow 
water.  It  appeared  to  have  the  essential  physics  capabilities  and  FE  techniques  for  this  kind  of  work, 
along  with  a  powerful  capability  permitting  users  to  incorporate  their  own  mathematical  expressions, 
of  almost  any  complexity,  in  almost  any  part  of  the  code.  It  is  this  latter  feature  that  has  made  it 
possible  to  add  the  results  of  the  above  R&D  to  the  Comsol  software.  Now,  at  the  end  of  FY08,  it  is 
clear  that  the  code  has  fulfilled  its  promise,  perfonning  well  beyond  initial  expectations.  It  is  an 
excellent  platfonn  for  both  R&D  and  production  modeling,  making  possible  all  the  research  advances 
and  modeling  results  described  above.  Looking  ahead  to  FY09  and  beyond,  and  considering  R&D 
currently  underway,  it  seems  reasonable  to  expect  that  the  code  will  be  able  to  handle  essentially  all 
the  modeling  required  on  these  ONR  programs. 

RELATED  PROJECTS 

SERDP  (Strategic  Environmental  Research  and  Development  Program). 
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Assessing  Sonar  Performance  Against  Underwater  UXO.  Apply  littoral  target  scattering  modeling 
tools  described  herein  to  underwater  UXO,  focusing  on  characterization  and  remediation  of  UXO- 
contaminated  environments. 
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